Coherent forward scattering peak induced by Anderson localization 
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Numerical simulations show that, at the onset of Anderson localization, the momentum distribu- 
tion of a coherent wave packet launched inside a random potential exhibits, in the forward direction, 
a novel interference peak that complements the coherent backscattering peak. An explanation of 
this phenomenon in terms of maximally crossed diagrams predicts that the signal emerges around 
the localization time and grows on the scale of the Heisenberg time associated with the localization 
volume. Together, coherent back and forward scattering provide evidence for the occurrence of 
Anderson localization. 
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Interference phenomena surviving statistical averages 
in spatially random media are of primary importance in 
order to understand the transport properties of bulk mat- 
ter. Already in weakly disordered media, spectacular 
deviations from the conventional Boltzmann picture of 
transport have been discovered and studied. Important 
examples include weak localization and universal con- 
ductance fluctuations in mesoscopic electronic systems, 
or coherent backscattering (CBS) and intensity correla- 
tions in speckle patterns in the context of wave transport 
in disordered media In strongly scattering systems, 
transport can even be completely suppressed, due to An- 
derson localization (AL) 0, H[. For several years now, 
this phase-coherent inhibition of transport by disorder 
has been under active scrutiny in widely different sys- 
tems using light Q, ultrasound [p | , microwaves [f|, and 
more recently ultracold atoms j7H12|. The observation 
of AL, however, remains a subtle issue, which requires 
accurate experiments, with a rather demanding control 
over absorption and dccoherencc. 

In order to closely monitor the localization dynamics, 
we have recently proposed to study the momentum dis- 
tribution of a coherent wave packet launched with finite 
velocity inside a random potential [l3|]. Phase-coherent 
multiple scattering is then clearly evidenced by the CBS 
peak, originating from the constructive interference of 
counter-propagating scattering amplitudes [HI, [Hj]. At 
the onset of AL, the amplitudes cease to propagate, so 
that the CBS peak should become stationary, an effect 
that we suggested to use as an indicator for AL [l3| . 
The issue proves in fact much richer. In this Letter, we 
report numerical evidence that AL indeed leads to a sta- 
tionary CBS signal, but also triggers the appearance of 
a novel, coherent forward scattering (CFS) peak in the 
opposite, forward direction. As a net effect, the momen- 
tum distribution shows a remarkable twin-peak structure 
(see Fig. [1]). We propose a theoretical explanation of 
CFS via a combination of the maximally crossed dia- 




FIG. 1. (color online) Ensemble-averaged momentum dis- 
tribution n(kx,k y ,t) of a simulated matter wave launched 
with initial momentum feo = (f. 5/^,0) inside a 2D ran- 
dom potential with correlation length f. The time unit is 
tq = mC, 2 /h. (a) t = 10r^: multiple scattering generates a 
coherent backscattering peak at —feo- (b) t = 90^: a coher- 
ent forward scattering (CFS) peak has emerged at +fco. (c) 
t — 950-Tf: the CBS and CFS twin peaks in the Anderson 
localization regime. 



grams responsible for CBS. This theory implies that the 
CFS peak, while almost indiscernible for weak disorder 
in the diffusive regime, starts to rise around the localiza- 
tion time Tioc = £ 2 /D (£ is the localization length and 
D the diffusion constant determined from the initially 
measured transport mean free path and time). Further- 
more, we find that the CFS contrast grows on the scale 
of the Heisenberg time th = hv^ d associated with the 
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FIG. 2. (Color online) Contrast of CBS (red circles) and CFS 
peak (blue squares) as a function of time, as obtained by 
fitting the data of Fig. [T]with the theoretical shapes, Eqs. ([2]) 
and 0, and a running average over a time window of 10t( . 



localization volume, where h is Planck's constant and v 
the average density of states per unit volume. The CBS 
and CFS twin peaks are thus a clear signal for the occur- 
rence of AL, and studying their dependence on system 
parameters and external fields promises further insights 
into wave dynamics in strongly disordered materials. 

We consider, then, a quasi-monochromatic wave packet 
that is prepared at time t = with mean wave vector 
fco and small spread Ak <C |fco| inside bulk disorder. 
In current cold-atom experiments fliil [l7j , such a wave 
packet is produced with ultracold atoms released from a 
trap inside an optical speckle field [lU ■ The atoms can be 
suspended against gravity by a magnetic field gradient, 
and accelerated to the desired fco by & suitable magnetic 
kick. After some time of evolution t inside the disorder, 
all fields are switched off, and the resulting momentum 
distribution n(fc,i) is accessed by time-of-fiight imaging. 

We first study the matter wave dynamics by numeri- 
cally solving the Schrodingcr equation with Hamiltonian 
H = p 2 /(2m) + V(r) in two dimensions (2D). Without 
loss of generality, V(r) — 0, where the overbar denotes 
the ensemble average over disorder realizations. The pair 
correlation function V(r)V(0) ~ V 2 C(r/() defines the 
potential strength via the variance V 2 . It also speci- 
fies the spatial correlation length which defines a time 
scale T£ = m( 2 /h and energy scale Eq = ti 2 /{mC, 2 ) [l9|. 
Fig. [T] shows the momentum distribution at three dif- 
ferent times, for a rather strong potential V = 5E^, 
averaged over 960 realizations of the speckle potential. 
The initial conditions are fco = (1.5/£, 0) and Ak = 
\/2fco/60 « 0.0236fco. At short times, in the diffusive 
regime [Fig. [lja) for t = 10tq], the CBS peak at —fco is 
the only distinctive feature. At longer times [Fig. [ljb) 
for t = 90t^], a CFS signal has risen at +fco- The CBS 
and CFS peaks then evolve toward a profile that remains 
stationary on the longest times accessible numerically 
[Fig. [TJc) for t ~ 950tq\. The peak contrasts, defined 
as height over background, are shown in Fig. [5] as a func- 
tion of time. The rise of the CFS signal goes along with a 
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FIG. 3. Constructive interference of amplitudes along 
maximally-crossed trajectories, leading to peaks at (a) —fco 
for CBS and (b) +fco for CFS, respectively. 



reduction of the CBS contrast. In the long run, both con- 
trasts converge to the same value. Similarly, the widths 
become comparable, and at the end of the time evolu- 
tion covered here, the coherent signals appear as twin 
peaks. This progressive symmetrization of the momen- 
tum distribution is consistent with the expectation that 
it should become invariant under fc i— > — k at very long 
times because of time-reversal symmetry [2Cj| . 

The remarkable symmetry between the two peaks 
suggests that CFS has the same physical origin as 
CBS, namely the interference of amplitudes counter- 
propagating along the same scattering path [2l|. This 
interference contribution is a soft mode of central impor- 
tance, known as the cooperon in mesoscopic systems 
Fig. OUa) represents the real-space trajectories associated 
with CBS. As shown in Fig. [3Jb), a CFS signal can ap- 
pear by concatenating two CBS processes, flipping twice 
the initial momentum and leading to an enhancement 
around +fco- Higher-order CFS (CBS) contributions can 
arise from chaining an even (odd) number of maximally- 
crossed diagrams. The remaining question then is to un- 
derstand why the CFS signal, while small in the diffusive 
regime, is boosted in the Anderson localization regime. 

The quantitative analysis is simplified by assuming a 
very narrow initial momentum distribution n(fc, t = 0) = 
(27r) rf (5(fc — fco) (see below for the impact of a finite mo- 
mentum spread Afc > 0). The ensemble-averaged mo- 
mentum distribution at a later time t, 



n(k, t) = 



dE 



dw 
27r' 



(1) 



is determined by the kernel <&kk E of intensity propaga- 
tion from initial momentum fco to final momentum fc via 
energy E. This kernel admits a well-known perturba- 
tion expansion in the disorder regime The isotropic 
diffusive background is described by the series of lad- 
der diagrams represented in Fig.^a), which results in a 
maximum background value of T s j{hmv) = riL, where r s 
is the elastic scattering time (l3j . 

The CBS contribution on top of the background is gen- 
erated by the series of maximally-crossed diagrams of 
Fig.[3ja). Its central ingredient is the diffusive cooperon 
Pb(w, q) = [—iuj + D(E,uj)q' 2 }^ 1 that peaks around 
q = fco + fc = 0. In the diffusive regime t t\ oc , D(E) is 
independent of uj. The resulting CBS signal then has a 
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FIG. 4. Ladder (a) and maximally crossed (b) multiple scat- 
tering Feynman diagrams, giving rise to the diffusive back- 
ground and the CBS signal, Eq. ([2]). Diagram (c) is respon- 
sible for the CFS peak, Eq. (J5]). Diagram (d) is exactly zero. 



decreasing width (2Dt) -1 / 2 , because the interfering am- 
plitudes originate from points whose distance grows dif- 
fusively [lH]. In the localization regime t ^> ri oc , this dif- 
fusive expansion stops, and the diffusion constant obeys 
the scaling relation D(E,w) ~ —iuj^ 2 (E) (22[. Conse- 
quently, the CBS signal becomes a stationary Lorentzian 
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whose width is determined by the localization length 
£(Eo), evaluated at the effective energy Eq of the ex- 
citation with wave vector fco inside the random medium. 
The second term in Eq. ([2]) is a negative background 
contribution that ensures particle-number conservation, 
J dkn^(k) = 0. Technically, it is obtained by dressing 
the cooperon of Fig. [3ja) with additional impurity lines 
[l[ , commonly represented by the shaded 4-point Hikami 
box of Fig.lUV) [13. 

Let us now repeat this analysis for the CFS peak, which 
originates from the trajectories depicted in Fig.^b), in- 
volving the direct connection of two cooperons. Its con- 
tribution to the propagation kernel reads 
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J J^jd A (Q, E)P E (q + k^)P E (q + k ,uj), 
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(27T) C 

where A{k, E) is the spectral function and d the spa- 
tial dimension. In the diffusive regime, a straightforward 
calculation shows that the CFS peak has a contrast in 
2D of the order of l/(ko£), where I denotes the transport 
mean free path. Thus, CFS is too weak to be visible in 
the weak-disorder regime k$£ 3> 1. In the localization 
regime, however, the scaling D(oS) ~ —icu^ 2 makes both 
cooperons singular and results in a signal that grows lin- 
early in time: 



n( cc )(fc,i) = 3n L — F d (\k-k \0, 

TH 



(4) 



with Fi (x) 



1/(4 



x 2 ), F 2 (x) 

asmh(x/2)/(nxy/4 + x 2 ), and F 3 (x) = atan(x/2)/(47ra;) 
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FIG. 5. Early-time CFS contrast for disorder strengths V = 
5E( (red circles), V = 4£^ (green squares), V = ?>E^ (blue 
triangles), and V = 2E^ (purple diamonds). According to our 
theory, the CFS begins to rise when localization sets in. The 
inset shows ti oc , estimated by the earliest CFS appearance, 
as a function of V. 



This CFS signal is indeed an interference peak centered 
at +fco whose width is again given by the localization 
length, and whose contrast is of the order of t/rji- 

Just as the CBS signal, the CFS contribution should 
also conserve the particle number. This leads us to dress 
the double cooperon with the 6-point Hikami box shown 
in Fig. HJc). This diagram also appears in second-order 
WL corrections to the electronic conductivity (23|. As a 
result of summing 16 different contributions, Eq. (|4]) is 
modified to: 



r(CC), 



t 



n^~>(k,t) = n L —F d {\k - feolO +n^>(k,t), (5) 

TH 

where n^ 2 \k,t) is a small negative background con- 
tribution that guarantees particle-number conservation, 
J dfcn' cc )(fc, t) = 0. In principle, there could be other 
isotropic contributions with an order of magnitude com- 
parable to the CFS signal, such as the one depicted in 
Fig.H{d), a combination of ladder and crossed series. But 
this contribution is exactly zero after proper dressing by 
the Hikami box. Similarly, other possible terms of the 
same order of magnitude vanish. This clearly identifies 
the class of diagrams |4jc) as the one responsible for the 
CFS peak in the Anderson localization regime. 

We briefly evaluate the effect of a finite source coher- 
ence: the convolution of the signal (|5|) by a source dis- 
tribution with width Afc > reduces the contrast by a 
term of order (Afc£) 2 . From our numerics for V = 5E^, 
we estimate £ ~ 8( from the width of the CBS peak, such 
that the reduction remains small, (Afc£) 2 « 0.1. Having 
a coherent enough source with Akt; 1 is thus manda- 
tory to observe the CFS peak. In the opposite case, the 
CBSpeak will have already decayed by the time AL sets 
in [13j, now without any remarkable signature on the 
isotropic background. 

Eq. J5]) implies that from t ~ t\ oc onwards, the CFS 
peak rises to a contrast of order unity on the Heiscnbcrg 
time scale th = hv£ d associated with the localization 
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volume. Fig. [5] shows the growth of CFS at early times 
obtained numerically for different values of V. Using the 
CFS onset as a measure for the localization time, we show 
in the inset of Fig. [5]ti oc as a function of V. The observed 
decrease agrees with the expectation that the localization 
time is smaller for stronger disorder. For the present pa- 
rameters, spectral broadening is very important (cf. the 
broad background in Fig. [T]). Consequently, the full sig- 
nal is the sum of contributions from a rather wide range of 
energies, with lowest-energy components being localized 
first. A fully quantitative theory for the strong-disorder 
regime would require to know spectral densities, diffu- 
sions constants and localization lengths, which are not 
yet at hand for spatially-correlated potentials. 

Our numerics shows that at long times t > th, the CBS 
and CFS contrasts saturate at a common value smaller 
than unity. This dynamics cannot be captured by our 
theoretical approach because it would require a summa- 
tion of the full CBS/CFS series of repeated cooperons to 
all orders, including the dressing by higher-order Hikami 
boxes, which poses a formidable task. The failure of 
perturbative expansions beyond the Heisenberg time tr 
is anyhow well known, and one has to resort to other 
tech niq ues like supersymmetry or random-matrix theory 
p4l [25| . We expect, however, that the two salient fea- 
tures of our explanation, namely (i) the emergence of 
CFS around the localization time t\ oc and (ii) a CBS 
and CFS dynamics on the Heisenberg time scale th, re- 
main valid. Indeed, a similar dynamics was predicted for 
the dynamical echo within weakly disordered quantum 
dots y26| and observed with ultrasound [27| ■ In our setup 
of bulk disorder, the confinement is due to AL, and the 
localization length stands for the effective system size. 
To our knowledge, it is an open question whether the 
CBS/CFS peaks can remain stationary or perhaps de- 
cay to a completely isotropic distribution on very long 



time scales associated with the glassy dynamics of the 
localization regime [28J. 

Future work should also address the impact of finite- 
size effects on the CFS signal. In many experiments, 
waves are sent through a scattering medium with open 
boundaries Especially for electromagnetic fields, 

quasi-plane- wave sources and detection of radiation pat- 
terns in the far field are readily available, enabling the 
clear observation of CBS in the reflected intensity [l| . We 
expect the CFS effect to be present as well in the trans- 
mitted intensity across a strongly localizing medium, but 
possibly difficult to observe because this requires good 
momentum resolution on the small transmitted signal. 

In conclusion, we have discovered a coherent interfer- 
ence phenomenon that is induced by Anderson localiza- 
tion. It manifests itself as an interference peak with high 
visibility in the forward direction and is closely connected 
to the CBS peak. We have proposed a theoretical expla- 
nation by a combination of maximally-crossed diagrams, 
and predicted CFS to occur in all dimensions d = 1,2,3. 
This effect provides an avenue for detailed investigations 
of AL in strongly random media, and should be in reach 
for current experiments with a good resolution in momen- 
tum space, such as cold atoms, excitons [29j . polaritons 
[30j , and photons HH . 
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